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The earlier introduced Differential HBT method aiming for detecting rotation is used to analyse 
the fluid dynamical model results of ultra-relativistic heavy ion reactions where the initial state 
has substantial angular momentum. The rotation effect and Kelvin Helmholtz Instability, lead to 
space-time momentum correlations, which can be detected by the method. 
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Peripheral heavy ion collisions have high angular mo- 
mentum, rapidly increasing with increasing beam energy. 
This is realized in rotating flow, large velocity shear, vor- 
ticity and circulation. Viscous, explosive expansion leads 
to the decrease of vorticity and circulation with time, 
however, with small viscosity the vorticity may remain 
significant at the final freeze-out (FO) stages. The ear- 
lier introduced Differential HBT method, a combination 
from standard two particle correlation functions, is ad- 
equate to analyze rotating systems. This was shown in 
few source models pQ. Now we look at this method and 
its results in a high resolution, Particle In Cell (PIC) 
fluid dynamics model. This model was used to predict 
the rotation in peripheral ultra-relativistic reactions [2], 
to point out the possibility of Kelvin Helmholtz Instabil- 
ity (KHI) [3], flow vorticity [I] and polarization [5]. The 
model was also tested for its numerical viscosity and the 
resulting entropy production [6]. 

One of the important observation in Ref. [ill is that in 
case of a globally symmetric fluid dynamical configura- 
tion the correlation function only includes cos(oi'u) and 
cosh(av) terms, therefore it will not depend on the direc- 
tion of the velocity, only on its magnitude. The direction 
dependence becomes apparent in the correlation function 
only if we take into account that due to the radial ex- 
pansion and the opacity of the strongly interacting QGP, 
the emission probability from the far side of the system 
is reduced compared to the side of the system facing the 
detector. 

Correlation Function from Fluid Dynamics: The 

pion correlation function is defined as the inclusive two- 
particle distribution divided by the product of the inclu- 
sive one-particle distributions, such that [TJ: 



C{pi,p 2 ) 



P2(j>l,P2) 
Pl{pi)Pl{p2Y 



(1) 



where p\ and P2 are the 4-momenta of pions. We intro- 
duce the center-of-mass momentum k = |(pi +P2) > 



1 The vector k is the wavenumber vector, k = p/h so for numerical 
calculations we have to use that he = 197.327 MeV fm., The same 
applies to q. 



and the relative momentum q = p\ — p 2 , where from the 
mass-shell condition [7J q Q = kq/k°. We use a method 
for moving sources presented in Ref. [8]. 



C(k,q) = l 



where 
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d X\d X2 cos[q(xi — x 2 )} x 



S(x 1 ,k + q/2)S(x 2 ,k-q/2) . 
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Using the emission function S(x, k), discussed in refs. [T], 
here R(k, q) can be calculated [5] via the function 



J(k,q) 



d 4 x S(x, k + q/2) exp(iqx) , 



(4) 



and we obtain the R(k, q) function as R(k, q) — 
Re[J(k,q)J(k,-q)}. 

We estimate the local pion density by the specific en- 
tropy, cr(x), as n w (x) cx n(x)a(x), where n(x) is the 
proper net baryon charge density. Thus the local invari- 
ant pion density is given by the Jiittner distribution as 



f J (x,P) 



n(x)o~(x) 



exp 



P^u^x) 
T(x) 



(5) 



where C„ = 47rm 2 r T '^(m^ /T), at temperature T, and 
K2 is a modified Bessel function. 

We asume that the single particle distributions, 
f(x,p), in the source functions are Jiittner distributions, 
which depend on the local velocity, u^{x), and we use 
the notation u\ = u(x\) = u^(x\). 

If we assume that the two coincident particles origi- 
nate from two points, x\ and x 2 , the expression of the 
correlation function, Eq. ^ will be become [T] 

R(k,q) = d 4 Xid 4 x 2 S(xi, k)S(x2, k) cos[q{x\— x?)] x 



exp 



u(xi) u(x 2 ) 
T( Xl ) ~ T{x 2 ) 



and the corresponding J(k, (j)-function is 

q ■ u(x) 



J(k,q) = / d x S(x,k) exp 



2T(x) 



(6) 

exp(iqx) , (7) 
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In Ref. Q] different one, two and four source systems 
were tested with and without rotation. Here we study 
only the case where the emission is asymmetric and dom- 
inated by the fluid elements facing the detector. 

In numerical fluid dynamical studies of symmetric 
(A+A) nuclear collision the initial state is symmetric 
around the center of mass (cm.) of the system, and (if 
we do not consider random fluctuations) this symmetry 
is preserved during the fluid dynamical evolution. 

Let us consider the usual conventions, z is the beam 
axis, and the positive z-direction is the direction of the 
projectile beam. The impact parameter vector points 
into the positive ir-direction, i.e. towards the projectile. 
Finally the y-axis is orthogonal to both. 

The fluid dynamical system, without fluctuations can 
be considered as a set of symmetric pairs of fluid cells. 

The emission probabilities from the two fluid cells of a 
source pair are are not equal. 

If we have several fluid cell sources, s, with Gaussian 
space and time emission profiles, then the source function 
in Jiittner approximation is 



2N3/2 7s" S (z) (^s) 



(2nR 



E 



c s 



exp 



(8) 



where n s =71^, and the spatial integral over a cell volume 
is, V ce ii = (2nR 2 ) 3 / 2 while the time integral is normalized 
to unity. Similarly the J-function is 



J(k,q) 



d A x S s (x,k)e iqx . (9) 



We then assume that the freeze-out (FO) layer is rela- 
tively narrow compared to the spatial spread of the fluid 
cells, so that the peak emission times, t s , of all fluid cells 
are the same. Then the exp(io^ s ) factor drops out from 
the J(/c, q)J(k, —q) product, rj This FO simplification 
is justified for rapid and simultaneous hadronization and 
FO from the plasma. For dilute and transparent matter 
the correlations from the time dependence of FO should 
be handled the same way as the spatial dependence. 

Due to mirror symmetry with respect to the [x, z], re- 
action plane, it is sufficient to describe the cells on the 
positive side of the j/-axis. The other side is the mirror 
image of the positive side. Then we can evaluate the 
correlation function the same way as in Ref. pQ. 



2 If the emission is happening through a layer with time-like nor- 
mal, but the peak is not at constant t s , but rather at constant 
t s , then we can adapt the coordinate system accordingly, i.e. we 
can use the r, r] coordinates instead of t, z, see e.g. [9]. 



Thus we define the quantities: 



= (2^i? 2 ) 3/2 exp 
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where u® = j s , the local 4-direction normal of the mean 
particle emission from an ST point of the flow is (as- 
sumed to be time-like), R is the size (radius) of the fluid 
cells, and O s is the path length of the time integral from 
the ST point of the source, s, while assuming a Gaussian 
emission time profile pQ. The weights, w s arise directly 
from the Cooper- Frye formula [9]. 

We can reassign the summation for pairs, so that s = 
{i,j, k} will correspond to a pair of cells: at {i, j, k} and 
its reflected pair across the cm. point at the same time 
at {i*,j*,k*y. Then the function S(k,q) becomes 



d 4 xS(x,k) = (2nR 2 f /2 x 
few. 



w s exp 



w s exp 



T s ) \ T, 

while, the function J{k, q) becomes 



(11) 



J(k,q) = Q c Y J Ps 



Qy'w s exp 
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+QMw* s exp 



(12) 



Only the mirror symmetry across the participant cm. 
is assumed, which is always true for globally symmetric, 
A+A, heavy ion collisions in a non-fluctuating fluid dy- 
namical model calculation. Then the correlation function 
can be evaluated using Eqs. ([2j(4]) . 

Based on the few source model results the Differential 
HBT method was introduced by evaluating the difference 
of two correlation functions measured at two symmet- 
ric angles, forward and backward shifted in the reaction 
plane in the participant cm. frame by the same angle, 
i.e. at t] = iconst., so that 



AC(k,q) = C(k + ,q out ) - C(k_,q out ). 



(13) 



For the exactly ±x -symmetric spatial configurations (i.e. 
k +x = k_ x and k +z = —k_z), e.g. central collisions or 
spherical expansion, AC(fc, q) would vanish! It would 
become finite if the rotation introduces an asymmetry. 

The Freeze-out Weights: While the flow symme- 
tries discussed above depend on the initial collision sym- 
metries, the weights do not follow these. The ansatz 
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used in Ref. [T] leads the weights as given in Eq. (10). 
These weights depend on the local mean emission direc- 
tion <t m , the flow velocity at the emission point and the 
opacity along the integrated path of the propagation of 
the emitted particle along this direction. The emission 
probability from different points of the space time is dis- 
cussed in more detail in Ref. jL based primarily on Refs. 
[lOUllj . These considerations may be more involved than 
the ansatz taken over from |5] . The determination of the 
FO surface normal or the mean emission direction from 
the ST freeze-out layer and the emission profile in this 
layer are the subjects of present theoretical research, see 
[T21 - fT5] . Here we do not discuss the effect of opacity along 
the path of the emission. 

The detector configuration is given by the two particles 
reaching a given detector in the direction of k. Thus the 
emission weights depend on the direction of the normal 
of the emission surface and of the emission, i.e. a and k. 
Furthermore, most monte-carlo cascade type of studies 
indicate [HI E3 HI] , that the majority of particles freezes 
out in a layer along a constant proper time hyperbola, 
with a dominant flow 4-velocity, which is normal to this 
hyperbola: <7 M ~ The origin of the hyperbola is 

at a ST point, which at low beam energies precedes the 
impact of the Lorentz contracted nuclei [Hj . 

We assume in the actual numerical calculations that in 
the expression of the weight, Eq. ( fl0| , Q s (q) is the same 

for all surface layer elements: — and 6 S = 9, so 
that w s = (k„ exp(-6 2 g 2 /2) , where & sfl = (cr°,er s ), 
so that = k Q a® + kcr s . If emission path time- 

length, 0, tends to zero, then the time modifying factor 
becomes unity. With the choice <7„ = ?i„, the time-like 



FO normal is cr s 



(7»,««). Then (k^) 



ry s k +ku s 



So the weight becomes 

(•j s k +ku s ) exp(-9 2 <7o/2). 



(14) 



This weight is explicitly different for the mirror im- 
age cell at x* — » — x s , where u* s — > —u s and then 
w* s = (7 s fc -fcu s )exp(-9 2 g 2 /2) . 

The weight factors appear both in the nominator and 
denominator of the correlator, so its normalization is bal- 
anced. On the other hand the role of the different factors 
in the weight have an effect to determine which cells con- 
tribute more or less to the result. This is a fundamental 
problem of the FO process, and in real situations it can 
get further convoluted if the hadronization and FO coin- 
cide in a rapid expansion. 

Results: The sensitivity of the standard correlation 
function on the fluid cell velocities decreases with de- 
creasing distances among the cells. So, with a large num- 
ber of densely places fluid cells where all fluid cells con- 
tribute equally to the correlation function, the sensitivity 
on the flow velocity becomes negligibly weak. 

Thus, the emission probability from different ST re- 
gions of the system is essential in the evaluation. This 
emission asymmetry due to the local flow velocity occurs 
also when the FO surface or layer is isochronous or if it 
happens at constant proper time. 
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FIG. 1. (color online) The dependence of the standard cor- 
relation function in the fc+ direction from the collective flow, 
at the final time. 



We studied the fluid dynamical patterns of the calcula- 
tions published in Ref. [3] , where the appearance of the 
KHI is discussed under different conditions. We chose 
the configuration, where both the rotation [2J, and the 
KHI occurred, at b = 0-7b max with high cell resolution 
and low numerical viscosity at LHC energies, where the 
angular momentum is large, L m 10 e h [13] . 

We used the Differential HBT method introduced in 
Ref. [T] for simplified examples of a few fluid cells. 
These examples were spatially symmetric, thus the stan- 
dard correlation function did not show any difference if 
it is measured at two symmetric k and q-out angles, e.g. 
in the reaction, [x-z] plane at k + — (k x ,0, +k z ) q + = 
(q x ,0,+q z ) and k = (k x ,0,-k z ) q = (q x ,0,+q_), i.e. 
AC(k,q) vanished. Here we have chosen two directions 
at rj = ±0.76, that is at polar angles of 90 ± 40 degrees. 
These are measurable with the ALICE TPC detector and 
at the ATLAS and CMS detectors also. 
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FIG. 2. (color online) The differential correlation function 
AC(k, q) at the initial and final times of flow development. 

The standard correlation function is both influenced 
by the ST shape of the emitting source as well as its 
velocity distribution. The correlation function becomes 
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narrower in q with increasing time primarily due to the 
rapid expansion of the system. At the initial configura- 
tion the increase of |fe| leads to a small increase of the 
width of the correlation function 
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FIG. 3. (color online) The flow velocity dependence of the 
differential correlation function at the final time. 

Nevertheless, in theoretical models we can switch of 
the flow, and analyse how the flow influences the corre- 
lation function and especially the differential correlation 
function, AC(k,q). 

Fig. [T] compares the standard correlation functions 
with and without the collective flow at the final time 
moment. Here we see that the flow causes a decrease of 
the width in q for the distribution at high values of |fe|. 

The differential correlation function is more sensitive 
on the flow as it is shown in the few source models. In 
Fig. [2] AC(k, q) is shown for the initial and final times. 
Although the amplitude of the differential correlation 
function is smaller the differences between the initial and 



final configurations are significant. The dependence on 
|fe| is especially large at the final time. 

At the final time the flow dependence is again minimal 
at small momenta, |fe|, while at large momenta the flow 
leads to a change in the shape of the differential correla- 
tion function, and leads to an increase of its amplitude, 
Fig. [3] Here the expansion and the flow interfere with 
each other and this leads to a more complex behaviour. 
(The spatial configuration is not symmetric for ±x, and 
this leads to nonvanishing AC(fc, q) at v — in Fig.[3j) 

These tests demonstrate that the differential correla- 
tion function is sensitive to the flow and its rotation. 

Conclusion: The present initial studies verify the de- 
pendence of the differential correlation function on the 
flow in general. At the same time it also shows that in 
realistic model situations the effect of the spatial size and 
shape, the flow expansion and the flow rotation all influ- 
ence the differential correlation function. To entangle the 
three different effects from one another further more de- 
tailed model studies are necessary and a more detailed 
evaluation of the differential correlation function would 
be beneficial, e.g. in addition the present q — out studies, 
the q — side and q — long components of the correlations 
functions should also be studied. Furthermore studies 
where k points in the out of plane (y) direction may also 
be beneficial. 

In addition theoretical model studies should be com- 
plemented with experimental investigations. In this case 
it is important to determine the precise Event by Event 
cm. position of the participants |17| . to be able to mea- 
sure accurately the emission angles, which are crucial in 
the present AC(k,q) studies. 
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